Joint spectrum of photon pairs measured by coincidence Fourier spectroscopy 
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We propose and demonstrate a method for measuring the joint spectrum of 
photon pairs via Fourier spectroscopy. The biphoton spectral intensity is 
computed from a two-dimensional interferogram of coincidence counts. The 
method has been implemented for a type-I downconversion source using a pair 
of common-path Mach-Zender interferometers based on Soleil compensators. 
The experimental results agree well with calculated frequency correlations that 
take into account the effects of coupling into single-mode fibers. The Fourier 
method is advantageous over scanning spectrometry when detectors exhibit 
high dark count rates leading to dominant additive noise. 
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Correlated pairs of photons are a popular choice in 
efforts to implement emerging quantum-enhanced tech- 
nologies. Proof-of-principle experiments have demon- 
strated ideas such as quantum cryptographyji quan- 
tum clock synchronization, 2,3 quantum optical coherence 
tomography^ and one-way quantum computingA In par- 
al·lel with the expanding range of potential applications, 
the need to develop appropriate tools to engineer and 
to characterize sources of photon pairs is becoming ap- 
parent. Among various degrees of freedom describing 
optical radiation, the spectral one is essential to a num- 
ber of techniques. 2,3,4 Also in other protocols, based on 
degrees of freedom such as polarization^ the spectral 
characteristics needs to be carefully managed in order 
to ensure the required multiphoton interference effects. 
This demand has brought a number of methods to con- 
trol the spectral properties of photon pairs by engineer- 
ing nonlinear media and the pumping and the collection 
arrangements£ïi2à2à±íi A development that is needed to 
match these advances is the ability to diagnose accurately 
two-photon sources and to measure reliably their char- 
acteristics. An important work in this context is the 
recent application of scanning spectrometers to obtain 
joint spectra of photon pairsü 

In this paper we demonstrate experimentally two- 
photon Fourier spectroscopy as a method to measure the 
joint spectrum of photon pairs. The setup is based on two 
independently controlled Fourier spectrometers in the 
common-path configuration. Such an arrangement guar- 
antees long temporal stability necessary to characterize 
weak sources of radiation operating at single-photon lev- 
els. We show that the two-dimensional map of coinci- 
dence counts recorded as a function of delays in two in- 
terferometers can be used to reconstruct the joint spec- 
trum of photon pairs. We present a measurement for 
a type-I spontaneous down-conversion process in a bulk 
beta-barium borate (BBO) crystal, and compare the re- 
sults of the reconstruction with theoretical predictions. 



An idealized scheme of the experimental method is pre- 
sented in Fig. G] The source X produces nondegenerate 
pairs of photons in distinct spatial modes represented by 
diverging lines. Each photon is sent into a separate inter- 
ferometer, where it is divided on a beam splitter BS into 
two wavepackets subjected to a delay difference ta and 
tb for the photon A and B. The two wavepackets inter- 
fere on a beamsplitter BS and light emerging from the 
output of the interferometer is detected. A coincidence 
event is recorded if both photons reach their respective 
detectors. The quantity of interest is the coincidence 
probability pabÍJa-, Tb) measured as a function of the 
delays ta and tb- 

The measurement is carried out in the regime when the 
response time of the detector is much longer than the 
inverse of smallest bandwidth characterizing the spec- 
trum of the source. Then the only relevant character- 
istics of the source is the joint spectral intensity given 
by (:Ía(üüa)Íb(u b )··)i where Ïí(ují) (i = A,B) are the 
spectral intensity operators for the beams A and B, 
and (: ... :) denotes the quantum mechanical expecta- 
tion value of normally ordered operators. For a pho- 
ton with a well defined frequency üü, the probability 
of reaching the detector is given by the Standard ex- 
pression (1 + coscc;tí)/2. Consequently, the probability 
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Fig. 1. A simplified scheme of the double Fourier spec- 
trometer applied to a source X of photon pairs. The pho- 
tons enter interferometers with independently adjusted 
optical delays ta and tb and are counted by detectors 
and D#. 
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VabÍJa, tb) of a coincidence event for a source with an 
arbitrary spectrum takes the form: 



(a) 



Pab(t a ,t b ) cx - J dcj A J dcj B ('.Ia(wa)Ib(wb).) 

x(l + cosüúata)(1 + cosüj b tb) (1) 

An example of the coincidence interferogram is shown in 
Fig. Efa). The two-dimensional Fourier transform of the 
coincidence probability PabÍJa^ Tb) comprises the follow- 
ing terms: 

J dT A dT B PAB(TA, T B ) exp(iO A T^ + Ítt B T B ) 

oc 6(íl A )6(íl B )(:N A N B :) 

+ ^6{íl A ){: N A Í B (\íl B \) :) + ^S(íl B )(: Í A (\Ü A \)N B :) 
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(:Í A (\n A \)ï B (\Q B \):) 
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where Ni = J duo Ïí{oj) is the operator of the total pho- 
ton flux in the zth beam, i = A,B. The first term, 
localized at Q A = Ü B = 0, is proportional to the to- 
tal number of photon pairs. The two middle terms lie 
on the axes Ü A = or Ü B = 0, and have the shape 
of the single-photon spectra conditioned upon the de- 
tection of the conjugate photon. Finally, the last term 
contains the sought joint two-photon spectrum. For op- 
tical fields, these terms occupy distinct regions in the 
Üa^b plane and can be easily distinguished, as shown 
in Fig. Etb). It is helpful to trace the origin of the 
four terms on the right hand side of Eq. (0) to the co- 
incidence interferogram. The vertical and the horizon- 
tal fringes generate the terms 5(Q A )(: N a Ib{\^b\) : ) and 
S(Qb)(· Ia(\Qa\)N b :), whereas it is the diagonal fringe 
pattern that contains information about the joint spec- 
trum (: Ia(\^a\)Ïb(\^b\) '-)• This defines the region of 
the coincidence interferogram that needs to be scanned in 
order to compute the joint spectrum. It has the shape of a 
tilted rectangle outlined in Fig. Efa). Let us note that the 
grid spacing in a given direction can be adjusted to the 
characteristic scale of interferogram structures. Specifi- 
cally, the grid can be sparse in the direction parallel to 
the fringes, while in the perpendicular direction it needs 
to be fine enough to resolve the oscillations. Then the 
Fourier transform of the experimental data covers the 
region marked with a dashed rectangle in Fig.Efb), con- 
taining the joint spectrum. 

Our experimental setup is depicted in the Fig. CJa). 
The photon pairs were generated in a 1 mm thick non- 
linear BBO crystal in a type-I process. The crystal was 
pumped by 100 fs long pulses centered at 390 nm, 20 mW 
average power and a repetition of 80 MHz. The ultra- 
violet beam was focused on the crystal to a spot mea- 
sured to be 155 /im in diameter. The crystal was cut at 
29.7° to the òptic axis, and oriented perpendicular to the 
pump beam. Two Mach-Zender interferometers MZ1 and 
MZ2 collected down-converted light at angles 1.28° and 
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Fig. 2. (a) A typical coincidence interferogram as a func- 
tion of optical path differences ta and tb and (b) its 
Fourier transform. The dashed rectangle in (a) defines 
the scan range, which with suitable sampling density 
yields the region of interest in the frequency domain, out- 
lined in (b). 



1.05°. The photons transmitted through the interferome- 
ters were coupled into single-mode fibers defining the spa- 
tial modes in which the down-conversion is collectedi^ 
Finally the photons were detected using single photon 
counting modules SPCM connected to fast coincidence 
electrònics and a PC controlled counter board. 

In order to ensure the stability of the interferomet- 
ric setup over the entire two-dimensional scan we used 
a pair of common path Mach-Zenders, in which the two 
arms were implemented as orthogonal polarization com- 
ponents while the optical path difference was modulated 
by a Soleil compensator, as depicted in Fig. CJb). The 
generated photons entered from the left, with their polar- 
ization set to 45° by the half-wave plate HWP. Then the 
photons went through a block of crystalline quartz QP 
with a vertical òptic axis, and a pair of wedges QW with 
horizontal òptic axes. With this setup, the sign and the 
value of the delay between the horizontal and vertical po- 
larization components could be accurately controlled by 
sliding one of the wedges, mounted on the stepper motor 
driven translation stage, into the beam. Finally the two 
polarizations were brought to interference using a second 
half-wave plate HWP and a polarizing beam-splitter P, 
and the transmitted photons were coupled into a single- 
mode fiber F using aspheric Iens L and detected. We ver- 
ified that for the spectral range of interest birefringence 
dispersion could be neglected and consequently the de- 
lays tí were linearly dependent on the displacements of 
the quartz wedges. The quartz blocks used in the setup 
allowed to vary the delays from —250 fs to 250 fs. 

The complete measurement consisted in a scan of a 
rectangular grid depicted in Fig. [Ja) spanned by 800x40 
points, where the first number refers to the direction per- 
pendicular to the fringes. The corresponding mesh was 
0.57 fs x 2 fs, and coincidences were counted for 3 s at 
each point. The reconstructed joint spectrum of pho- 
ton pairs, normalized to unity, is depicted in Fig. E{a). 
We compare it with theoretical calculations plotted in 
Fig. Efb). The theoretical model used in these caleu- 
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Fig. 3. (a) Experimental setup. X, BBO crystal; MZ1, 
MZ2, Mach-Zender interferometers; SPCM, single pho- 
ton counting module. (b) Common-path Mach-Zender 
interferometer. HWP, half-wave plate; QP, quartz plate; 
QW, quartz wedges; P, polarizer; L, aspheric Iens; F, 
single-mode fiber. 
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a more quantitative comparison, Figs.^c) and (d) show 
the cross sections of the joint spectra along directions of 
maximum and minimum width in the frequency domain. 
In these plots, the experimental data have been inter- 
polated between the points of the Fourier-domain grid, 
and presented together with statistical errors calculated 
assuming Poissonian noise afïecting coincidence counts. 

In summary, we proposed Fourier spectroscopy for 
measuring the joint spectrum of photons pairs, and 
demonstrated its application to down-converted light 
generated in a type-I BBO crystal. The result of the 
reconstruction agrees well with a careful theoretical cal- 
culation of the joint spectrum. We were able to reduce 
substantially the overall duration of the measurement 
by selecting the region of the interferogram which con- 
tains information about the relevant characteristics of the 
spectrum. Finally, let us note that compared to scan- 
ning spectrometers, Fourier spectroscopy exhibits higher 
signal-to-noise ratio when detection noise is dominated 
by an additive contribution. This effect, known as the 
multiplex advantage, 13 is important in the case of high 
dark count rates, which is typical for single-photon mea- 
surement s performed at telecom wavelengthsi 
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Fig. 4. (a) Measured joint spectrum of photon pairs 
(' Ia(Xa)Ib(Xb) ') as a function of the wavelengths Xa 
and Xb] (b) theoretical predictions of the joint spectral 
intensity; the cross-sections of the theoretical (dashed 
line) and the experimental (sòlid with dots marking er- 
rors) distributions along the lines (c) A^ 1 + À^ 1 = const 



and (d) À y 



const passing through the maxi- 



lations assumed the exact phase matching function of 
the nonlinear crystal. The transverse components of the 
wave vectors for the pump and down-converted beams 
were treated in the par axial approximation. The joint 
spectrum was calculated for coherent superpositions of 
plane-wave components of the down-conversion light that 
add up to localized spatial modes defined by the collect- 
ing òptics and single-mode fibers. In order to facilitate 
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